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Assembly line balancing is a well-known problem in operations research. It con- 
sists in finding an assignment of tasks to some arrangement of stations where the 
tasks are executed. Typical objective functions require to minimize the cycle time 
for a given number of stations or the number of stations for a given cycle time. 

Common to all variants of assembly Une balancing are the precedence constraints 
among the tasks. Different inequalities for modeling these constraints have been 
proposed in the literature and lead to different formulations for problems such as 
the simple assembly line balancing problem (SALBP). We theoretically compare 
these formulations in terms of the strength of their linear relaxations. Moreover, we 
propose two new families of inequalities and a stronger formulation of the station 
limits for the SALBP and show that one of the resulting formulations dominates all 
existing ones. 

We evaluate the efficiency of these formulations experimentally and show that the 
proposed formulation produces equal or better results than existing formulations. 
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1 Introduction 



Assembly lines are an efficient way of producing large quantities of standardized products (Boysen et al., 
2008). The products are successively assembled on some arrangement of stations, most com- 
monly a linear sequence. On each station workers execute the subset of tasks assigned to that 
station. Due to technical restrictions the tasks usually have to respect precedence constraints and 
therefore cannot be assigned arbitrarily to stations. 

The station time is the total time needed to execute the tasks assigned to a station and the 
cycle time of an assignment is the maximum station time. In problems of type 1 the goal is to 
minimize the number of used stations for a given cycle time, while in problems of type 2 the 
goal is to minimize the cycle time for a given number of stations. Other problem types include 
the feasibility problem (type F) of finding some solution for a given number of stations and a 
given cycle time, and the problem of type E of maximizing the line efficiency, i.e., minimizing 
the product of the number of stations and the cycle time. 

The most simple and well studied problem of this kind is the simple assembly line bal- 
ancing problem (SALBP) (Salveson, 1955). In this problem each task has a fixed execution 
time and the stations are arranged in a linear order. The four variants mentioned above are 
known as SALBP- 1, SALBP-2, SALBP-F, and SALBP-E. All these problem variants are NP- 
complete (Scholl and Becker, 2006). 

In this paper we theoretically and experimentally compare several formulations for the SALBP- 

1 and SALBP-2, and propose two new families of valid inequalities for expressing the prece- 
dence constraints. 

The remainder of this paper is organized as follows. In the next section we present a formal 
definition and a mixed-integer linear model of the SALBP- 1 and the SALBP-2. Next, we dis- 
cuss different models proposed in the literature. In Section 3 we theoretically compare these 
models and two new families of valid inequalities for expressing the precedence constraints. A 
computational study is presented in Section 5 and we conclude in Section 6. 

2 Mathematical models for the SALBP 

In this section we review mathematical models for SALBP- 1 and SALBP-2 that have been pro- 
posed in the literature. Let (A'^, <) be a partially ordered set of tasks. For tasks i,] € N, we write 
/ < i if / < i and / ^ j and / < j if / immediately precedes j. Each task / G has execution time 
tj. Let S be the set of stations. In the following we suppose that for the SALBP-1 an upper bound 
m on the number of stations is given (lA^^I is such an upper bound), and that S = [m]'. For the 
SALBP-2 the number of stations m is known. In this case we set S = [m]. 

There have been three types of models proposed in the literature. Bowman (1960) (in the 
revised formulation of White (1961)) proposed two formulations, one based on binary impulse 
variables, indicating the station to which a task is allocated, and another based on time variables, 
representing the starting time of the tasks. Scholl (1999) proposed a formulation using binary 
step variables, which indicate that a task is allocated to some station or a preceding station. 



We write [ra] = { 1 , . . . , m} 
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Since the formulation using time variables has been found inferior by Pastor et al. (2007) we 
focus in the following on formulations using impulse and step variables. 

2.1 Formulations using impulse variables 

An allocation of tasks to stations can be formulated using variables ;c„ S {0, 1} indicating that 
task / € A'^ has been allocated to station s GS. Any feasible allocation has to satisfy the occurrence 
constraints 

^X,; = l V/GA^, (1) 

seS 

which make sure that every task is assigned to a station, the precedence constraints 

Xtj< ^ Xsi V/, eN,i<j,yt e S, (2) 

seS\s<t 

and the nondivisibility constraints 

;c,,-e{0,l} ^ieN,s&S. (3) 

These constraints have been first proposed by Bowman (1960). Their above form is due to White 
(1961). 

For a given cycle time c and an upper bound m on the number of stations, the SALBP-1 can 
be formulated as 

(BWl-1) minimize ^sjs, (4) 

seS 

subject to ^ tiXsi < cys G S, (5) 

ieN 

Equations (l)-(3), (6) 
^.£{0,1} Vs€5. (7) 

For the SALBP-2 the cycle time c is variable and the number of stations m is fixed. It can be 
formulated as 

(BWl-2) minimize c, (8) 

subject to Yj ^'-^v' ys€S, (9) 

ieN 

Equations (l)-(3), (10) 

cGR. (11) 

These formulations are due to Baybars (1986), which also presents a survey of different for- 
mulations. In the following two subsections we present alternative formulations for the station 
limits and the precedence constraints. 
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2.1.1 Station limits 



For a given cycle time c and a given number of stations m, Patterson and Albracht (1975) ob- 
served that a task / € A'^ cannot be assigned to a station earlier than 



He) 

or to a station later than 



jeN\j<i 



Li{c,m) =m+\ 



(12) 



(13) 



Let A,(c,ni) = [£',(c),L,(c,m)] be the set of admissible stations for a task / € A'^. Then we can 
add to the formulation of the SALBP-1 the constraints 

Xsi = V/ G A^, A; (c, m) , ( 14) 

where m is an upper bound on the number of stations. Similarly, if c is an upper bound on the 
cycle time c we can add the constraints 

x,i = V/ G N,\/s ^Ai{c,m) (15) 

to the formulation of the SALBP-2. 

These inequalities have been first proposed by Patterson and Albracht (1975). They can be 
strengthened as proposed by Pastor and FeiTcr (2009). For the SALBP-1, by taking into account 
currently used stations, we obtain 

XL,(c,s)j<y. ysesyieN. (16) 

These inequalities are valid, since if a station s is unused, we can assume that all later stations 
are also unused and therefore the latest possible station is now Li{c,s — 1) = Li{c,s) — 1. 

For the SALBP-2 we can model the cycle time explicitly to obtain better bounds on the sta- 
tions. Let c be a lower bound on the cycle time and let C = [c,c] be the set of admissible cycle 
times. Then we can represent the cycle time by 

c = Y.tr„ (IV) 

tec 

Y.rr = \, (18) 

tec 

r,G{0,l} t£C. (19) 

This allows to add the following valid inequalities to the formulation of the SALBP-2: 

Xei<l- £ rr yieN,Ei{c)<e<Ei{c), (20) 

teC\e<Ei{t) 

xii<l- £ rt \fieN,Li{c,m) <l <Li{c,m). (21) 

teC\Li{t,m)<l 

These inequalities are easily seen to be valid. Suppose, for example, that Lrec|(;<£,(f) '"f = 1 in 
equation (20). Then station e comes before the earliest possible station for the current cycle 
time, and therefore x^, = 0. Similarly, if L(ec|L,((,m)</ = 1 in equation (21) station I is after the 
latest possible station for the current cycle time, and therefore xu = 0. 
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2.1.2 Precedence constraints 

Patterson and Albracht (1975) have proposed to formulate the precedence constraints as 

^T,^^"} yi,jeN,i<j. (22) 

seS seS 

Thangavelu and Shetty (1971) observe that an alternative formulation of the precedence con- 
straints is given by 

^(m -s+\){x,i-x,j) > V/,7 eN,i< j. (23) 

seS 

2.2 Formulations using step variables 

An alternative formulation can be obtained using step variables Xsi € {0, 1} indicating that task 
/ E A'^ has been allocated to a station t £ S such that t <s. When using step variables, we obtain 
occurrence constraints 

= 1 V/ G A^, (24) 

and precedence consti^aints 

Xsi > Xs j yij eN,i<j. (25) 

Additionally, the step variables have to satisfy the continuity constraints 

Xsi < Xs+i.i yi €^ N,s & S,s < m. (26) 

To simplify the formulation, we assume that xoi = for all / € A'^. Then, the SALBP-1 can be 
formulated as 



(SCl-1) minimize ^^ys, (27) 

seS 

subject to - yseS, (28) 

ieN 

Equations (24)-(26), (29) 

^.€{0,1} V^e^. (30) 

and the SALBP-2 can be formulated as 

(SCl-2) minimize c, (31) 

subject to '^ti{xsi—Xs-i,i) < c E 5, (32) 

ieN 

Equations (24)-(26), (33) 

c G M. (34) 
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2.2.1 Station limits 

The station limits proposed for the formulations using impulse variables be can also applied to 
the formulation using step variables. Using the above notation we can add to the formulation of 
SALBP-1 

xe.{c)-u = yieN (35) 

^L,(c,m),/ = 1 V/ G (36) 

and to the formulation of SALBP-2 

XE.(c)-i,i = \fieN (37) 

XL.{c,,n),i = ^ V/eA^. (38) 

Similarly, the strengthened station limits in the step model are 

^L,M-i,, > l-J. V/EA^ (39) 

for SALBP-1 and 

xei<l- £ r, yieN,Ei{c)<e<Ei{c), (40) 

teC\e<Ei{t) 

xi-ij > £ r, V/ G N,Li{c,m) < I < Li{c,m) (41) 

teC\Li{t,m)<l 

together with equations (17)-(19) for the SALBP-2. 



3 New inequalities for the SALBP 

3.1 Station limits 

We can strengthen the station limits proposed by Pastor and Ferrer (2009) as follows. In inequal- 
ity (16), when a task cannot be assigned to station Lj{c,s) it also cannot be assigned to any later 
station. This justifies 

E ^ui<ys yseS,ieN. (42) 

u£S\u>Li{c,s) 

Similarly, in equations (20) and (21) when a task cannot be assigned earlier than station e this 
holds also for stations preceding e, and when a task cannot be assigned later than station I, this 
holds also for stations following I. Therefore we can strengthen these inequalities to 

£ ^„,-<l- £ r, yieN,Ei{c)<e<Ei{c), (43) 

ueS\u<e teC\e<Ei{t) 

and 

£ < 1 - £ rt V/ e N,Li{c,m) < I < Li{c,m). (44) 

ueS\u>l teC\Li{t,m)<l 



6 



3.2 Precedence constraints 



We propose two new formulations of the precedence constraints: 



yij £N,i<j,t <s 
yij eN,i< j,keS. 



(45) 
(46) 



Inequalities (45) are valid, since for / < j task / G N cannot be executed on a later station than 
task j G A^. 

Proposition 1 Inequalities (46) are valid for BW 1-1 and BW 1-2. 

Proof. If T,ses\s>kXsi = the inequality is trivially satisfied. Otherwise, task / is executed on 
station k or later. But since / < j task j has to be executed also on station or a later station, i.e., 

llseS\s>kXsj = 1- ■ 

Proposition 2 Inequalities (46) imply inequalities (22). 
Proof. Summing over all ^ G 5 yields, for all /, j ^ N,i< j. 



keSseS\s>k keSseS\s>k 




Proposition 3 Inequalities (23) are equivalent to inequalities (22). 



Proof. 



Y^{m-S+\) {Xsi - Xs ;•) > 



seS 

Y^{m-s+ \)xsi >Y^{m-s+ \)xsj 



seS seS 




which, in the presence of equalities (1), are equivalent to 



52 ■^■^.V! < J^^^sj- 



seS seS 



Proposition 4 Inequalities (46) imply inequalities (2). 
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Proof. From (46), in the presence of equality (1), and considering that all Xsj > 0, we obtain, for 
all /, j ^ N, i< j, and for all k £ S, 

seS\s<k seS\s<k-l seS\s>k 

and again in the presence of equality (1) 

seS\s<k 

<=^1- E Xsi<^-Xk-hj 

seS\s<k-l 

.•ieS\s<k- 1 

Therefore, for all f = ^ — 1, /c € 5 

seS\s<t 

We are missing the last inequality 

j ^ -"-if ) 
seS 

but this is trivially obtained in the presence of equality (1), since Y,ses^si = 1- * 
Proposition 5 Inequalities (2) /m/jfy inequalities (45). 

Proof. In the presence of equahty (1) we can rewrite inequalities (2), for all /, j eN,i< j,t G S, 
as 



Xtj < 1 - E 



SI 

seS\s>t 



se5'|.s>f 

Therefore, for all /, j £ N,i< j, t < s, 

Xsi -\~Xtj ^ £ Xf;i -\-Xtj < 1 . 

seS\s>t 

■ 

Proposition 6 Inequalities (46) dominate inequalities (22). 

Proof. Let A/^ = {a,Z7} and a < and suppose m = 3. Then, xi^ = = 1 /2, xi/, = 3/4, and 
x^^ = 1 /4 satisfies (22), but not (46). ■ 
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Table 1: Different formulations of the SALBP-1. 



Precedence constraints 


Station limits 




Name 


Equations 


Res. 


Name 


Equations 


Res. 


PA 


(22) 




1 






BW 


(2) 


on^m 


2 


(14) 




TS 


(23) 


2 

on 


3 


(14),(16) 


mn 


RC 


(45) 


2 

on 


4 


(14),(42) 


Tnn 


RC 


(46) 


2 — 2 

on m 








Base model 










Suffix 


Equations 


Var. 


Res. 






-1 


1,3,4,5,7 


m{n + 1) 


rn + n 







Proposition 7 Inequalities (46) dominate inequalities (2). 

Proof. Let = {a,b} and a<b and suppose m = 3. Then, xia = xih = X2b = x^a = 1/2 satisfies 
(2), but not (46). ■ 

Proposition 8 Inequalities (22) and (2) are incomparable. 

Proof. The example given in the proof of Proposition 6 satisfies (22), but not (2). The example 
from Proposition 7 satisfies (2), but not (22). ■ 

Proposition 9 Inequalities (22) and (45) are incomparable. 

Proof. The example from Proposition 6 satisfies (22), but not (45). The example from Proposi- 
tion 7 satisfies (45), but not (22). ■ 

In summary, inequalities (46) dominate, to the best of our knowledge, all inequalities proposed 
in the literature for formulating the precedence constraints. Tables 1 and 2 give a summary 
of models with different precedence constraints and station limits for the SALBP-1 and the 
SALBP-2 and their number of variables and restrictions. We obtain a different model for each 
combination of the precedence constraints, the stations limits and the base model. Observe that 
equations (14) and (15) do not increase the number of restrictions, but reduce the number of 
variables. In both tables o denotes the order strength of the instance. The order strength of an 
instance is at most (2) . 

Figure 1 shows the relationships between these models. 

4 Comparison of formulations using impulse and step variables 

In this section we demonstrate that RC implies SC under the variable conversion 

u^S\u<s 
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Table 2: Different formulations of the SALBP-2. 



Precedence constraints Station limits 



Name Equations 


Res. 


Name 


Equations 


Res. 


PA (22) 


2 

on 


1 






BW (2) 


on^m 


2 


(15) 




TS (23) 


2 

on 


3 


(15),(17)-(21) 


mn 


RC (45) 


on^m 


4 


(15),(17)-(19),(43),(44) 


mn 


RC (46) 


2 2 

on m 








Base model 










Suffix Equations 


Var. 


Res. 






-2 1,3,8,9,11 


nm + 1 


m + n 







RC 



B 



BW 




PA 










RC 




TS 



B 



A « 



B 



A is better than B 

A and B are incomparable 

A and B are equivalent 



Figure 1: Relationships between SALBP formulations. The relationships are due to different 
precedence constraints and are valid for SALBP- 1 as well as SALBP-2 formulations 
with the same station limits. 
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Using this conversion, since Xsi > for all 5 G 5 and / G N, the continuity constraints of SCI are 
satisfied because 



weS|M<s hgS|m<4+1 



and by equality (1) we obtain 



ueS\u<\S\ ueS 



i.e., the occurrence constraints hold. We can write the precedence constraints (25) as 

u£S\u<s uES\u<s 

n65|u>s uES\u>s 
ueS|n>.s HGS|«>.V 

which shows that the precedence constraints in both formulations are equivalent, and we can 
write equation (28) as 



ieN \ueS\u<s iieS\u<s-l ' 



E5| 

ieN 

which shows that the cycle time constraints in both formulations are equivalent. 

The implication also applies when inequalities for station limits are present. By equation (14) 
we have x„' = for 5 < Ei(c) and for s > Li(c,m) and therefore equation (35) 



s<Li{c,m) sES 



and equation (36) 



s<Ei{c) 



hold. By equation (42) we have 



seS\s<Li(c,s) seS\s>Li{c,s) 

and by equations (43) and (44) we obtain 

xei= £ ^.,-<l- £ rt ^ieN,Ei{c)<e<Ei{c) 

s€S\s<e t€C\e<Ei{t) 

and 

xi-i,i= Y Xsi = l- Y ^^i^ yi eN,Li{c,m) <l <Li{c,m). 

seS\s<l-l seS\s>l teC\Li{t,m)<l 

This shows that the dynamic station limits in RC imply the dynamic station limits in SC. 
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5 Computational experiments 



We have empirically evaluated the performance of the formulations of the SALBP-1 and the 
SALBP-2 presented in Tables 1 and 2. We limited the comparison to the best known formula- 
tions PA and BW from the literature and the new formulation RC for the precedence constraints 
combined with all four sets of equations for the station limits, for a total of 12 SALBP-1 and 
SALBP-2 formulations. 

The formulations have been tested on 269 instances of the SALBP-1 and 302 instances of 
the SALBP-2. These instances are available online (ALB, 2011). Details about them can be 
found in Scholl (1993). Currently the optimal value is known for all except one of the SALBP-1 
instances and 14 of the SALBP-2 instances. In the evaluations below solutions for instances 
without a known optimum were never considered optimal. 

All experiments were performed on a PC with an Intel Core 17 CPU running at 2.8 GHz and 
12 GB of main memory. We used the solver CPLEX 12.3.0 with standard options, except for 
a MIP optimality gap of 10^^, running in deterministic mode with two threads for a maximum 
time of 600 seconds. All computation times reported are in seconds of real time. Following 
Pastor and Ferrer (2009) we use in our experiments lower bounds on the number of stations 
m = LMl := [LieA^^i/c] and on the cycle time c = LCI := max{max,e/vf,-, [Liew^i/'^l } for 
the SALBP-1 and SALBP-2 respectively (Scholl and Becker, 2006) and an upper bound on the 
number of stations of m = min{2m, lA^^I} and on the cycle time of c = 2c. 

Tables 3 and 4 show the results for the formulations of the SALBP-1 and the SALBP-2, 
respectively. For each tested model, the first column presents the relative deviation of the LP 
bound from the optimal value in percent. The next three columns present the number of instances 
for which the branch-and-cut solver of CPLEX found a provably optimal solution, the average 
solution time for these instances, and the average number of nodes explored in the branch-and- 
cut tree. The next two columns present the number of instances for which an optimal solution 
was found, but could not be proven to be optimal and the average number of nodes explored 
for these instances. The following two pairs of columns present the same information for the 
instances which terminated with a feasible, but not optimal solution, and for the instances for 
which no feasible solution was found. The last column presents the number of instances where 
the formulation obtained the best value found over all 12 formulations of the same problem type. 

In our evaluation, when comparing solution values of different formulations we inform the 
/7-values obtained by a conservative, non-parametric paired sign test, where half of the ties were 
assigned to each sample (Dixon and Mood, 1946). 

In both problem types, the relative deviations of the LP relaxation from the best known value 
improve with stronger station limits. The new station limits produce the best bounds, but im- 
prove only slightly over the station limits of Pastor and Ferrer (2009). For the SALBP-1 they 
lead to more provably optimal and more best solutions, but not for the SALBP-2. In both prob- 
lem types the solution value does not improve significantly compared to other forms of station 
limits (for p = 0.05). The relative deviations are the same under different formulations of the 
precedence constraints. Therefore, the influence of the precedence constraints on the results, 
discussed below, is due to different relaxations and cuts of subproblems during the execution of 
the branch-and-cut algorithm. 

Concerning the SALBP-1, formulations BW?i and RCn obtain better results than the corre- 
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Table 3: Comparison of formulations of the SALBP-1. Time and number of nodes are averages 
over the instances of each category. 











Optimal 








Sub-optimal 










Proven 


1 


Unproven 


Feasible 


Infeasible 




Model 


Dev. 


# 


Time 


Nodes 


# 


Nodes 


# 


Nodes 


# 


Nodes 


Best 


PAl 


5,408 


158 


32,5 


25027 


24 


667091 


47 


372095 


40 


24038 


205 


PA2 


5,171 


175 


24,8 


38739 


22 


495261 


35 


454219 


37 


60031 


218 


PA3 


4,349 


184 


34,5 


41886 


18 


993955 


31 


507752 


36 


4963 


215 


PA4 


4,348 


186 


35,2 


43814 


18 


1020934 


33 


561443 


32 


13527 


223 


BWl 


5,408 


174 


52,5 


66862 


20 


166021 


75 


109242 







221 


BW2 


5,171 


178 


28,0 


17152 


20 


280644 


71 


101728 







230 


BW3 


4,349 


189 


32,8 


34982 


16 


195908 


64 


108796 







232 


BW4 


4,348 


190 


22,5 


36286 


21 


443164 


58 


137030 







247 


RCl 


5,408 


184 


37,7 


32217 


16 


397677 


69 


94612 







234 


RC2 


5,171 


188 


28,0 


30878 


18 


217022 


63 


80017 







242 


RC3 


4,349 


189 


30,5 


30624 


15 


326869 


65 


168141 







242 


RC4 


4,348 


197 


35,6 


22313 


17 


284820 


55 


204217 







260 



Table 4: Comparison of formulations of the SALBP-2. Time and number of nodes are averages 
over the instances of each category. 





Optimal 


Sub-optimal 


Proven 


1 Unproven 


Feasible | Infeasible 



Model 


Dev. 


# 


Time 


Nodes 


# 


Nodes 


# 


Nodes 


# 


Nodes 


Best 


PAl 


2,100 


164 


70,4 


77770 


17 


258272 


121 


437489 







199 


PA2 


2,100 


170 


55,3 


56314 


12 


549199 


120 


330365 







196 


PA3 


0,692 


187 


57,4 


55120 


9 


657220 


86 


381573 


21 


1037 


206 


PA4 


0.691 


184 


28,3 


45916 


11 


1036212 


93 


529207 


14 


1151 


206 


BWl 


2,100 


188 


49,3 


59523 


17 


430656 


97 


254536 







219 


BW2 


2,100 


192 


50,1 


79907 


11 


347230 


99 


334954 







220 


BW3 


0,692 


186 


31,3 


49803 


10 


714271 


97 


359153 


9 


106 


210 


BW4 


0,691 


188 


24,9 


47618 


10 


334468 


97 


321281 


7 


85 


213 


RCl 


2,100 


189 


34,1 


54955 


14 


399572 


99 


424456 







234 


RC2 


2,100 


186 


32,4 


62847 


15 


224409 


101 


451419 







232 


RC3 


0,692 


200 


41,9 


86384 


10 


441196 


92 


308320 







225 


RC4 


0,691 


196 


47,0 


65991 


8 


402442 


98 


349380 







224 
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spending formulation FAn (with p < 0.04 for BW?i and p < 0.02 for RC?i, I <n <4). They are 
able to prove and find more optimal solutions, always find a feasible solution, and more often 
find the best result. The computation times are comparable, except the higher execution time 
of BWl. Formulations RCn prove and find more optimal solutions and more best values than 
BW?i, but the difference in solution value is not statistically significant (for p = 0.05). From a 
practical point of view formulation RC4 could be considered best, since it provides the largest 
number of best solutions in a time comparable to the other formulations. 

Concerning the SALBP-2, the formulations BW« find and prove more optimal solutions than 
formulations PA?i, but the improvement in solution value is not statistically significant. The 
reason is that BW?i obtains worse solutions on the 28 instances proposed by Scholl (1993), 
for which both formulations find only feasible solutions. Formulations RC?i find and prove 
more optimal solutions than BW« and PA«, except for RC2, and always find a feasible solution. 
They obtain better solutions values than PAn (with p < 0.001). The improvement over BW« 
is marginally significant (p < 0.04 for 1 < « < 3, and p = 0.06 for n = 4). Formulations RCn 
also find most best values. The computation times of formulations RCn and BW« are less 
than those of PA«, except for RC4. Computation times of RCn and BW« are similar, but a 
comparison of the computation times on the 142 instances that all formulations could prove 
optimal shows that RC?i tends to increase computation time. The difference in solution value 
among formulations RC?i is not statistically significant and we can observe a trade-off between 
proving more solutions optimal and finding better solution values, where formulation RC4 tends 
to find and prove less optimal solutions. From a practical viewpoint one may prefer formulation 
RCl which find the largest number of best solutions in a short time. Formulations BW2, RCl, 
and RC2 were able to prove that the previously open instance Wee-Mag 23 has optimal solution 
value 65. 

In a computational study Pastor et al. (2004) found no significant difference between formula- 
tions PA and BW, but observe that formulation BW leads to shorter solving times. Our results do 
confirm that only for the SALBP-2. The different results may come from the difference between 
the used solvers (CPLEX 8.0 and CPLEX 12.3). Pastor and Ferrer (2009) find that the dynamic 
station limits (PA3) increase the number of provably optimal solutions over formulation PA2, 
which is corroborated by our findings. The solution values, on the other hand, do not decrease 
significantly, and the dynamic station limits tend make it more difficult to find feasible solutions 
in the SALBP-2. 

6 Conclusions 

We have proposed two new formulations of the precedence constraints for the SALBP, and 
shown that one of them dominates, to the best of our knowledge, the precedence constraints 
proposed in the literature. They are applicable to other assembly line balancing problems with 
precedence constraints. We have further given a variable conversion under which the model 
for the SALBP using impulse variables implies the model using step variables when using the 
newly proposed precedence constraints and the strengthened station limits. Comparing the new 
and existing models, we have provided a complete classification of the relationships between the 
models currently existing in the literature. 
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An experimental evaluation corroborates the theoretical comparison. The newly proposed 
precedence constraints can improve upon the constraints of Patterson and Albracht (1975) and 
Bowman (1960), finding and proving more solutions optimal, and finding more best values. A 
conservative statistical test shows that the improvement of the solution value is significant, with 
one exception for SALBP-1, where the newly proposed constraints and those of Bowman (1960) 
have to be considered practically equivalent. 
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